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Abstract 



We address the role of 7-ray astronomy in the investigation of nonthermal 
processes in the large scale structure of the universe. Based on EGRET upper 
>^ ' limits on nearby galaxy clusters (GCs) we constrain the acceleration efficiency of 

CR electrons at intergalactic shocks to < 1% than the shock ram pressure. That 
^ ' implies a contribution to the cosmic 7-ray background from intergalactic shocks 

of order 25 % of the measured level. We model spatial and spectral properties 
^ ■ of nonthermal 7-ray emission due to shock accelerated cosmic-rays (CRs) in GCs 

^ ' and emphasize the importance of imaging capability of upcoming 7-ray facilities 

Qh' for a correct interpretation of the observational results. GC observations at this 

O ' photon energy will help us understand the origin of the radio emitting particles, 

^ • the possible level of CR pressure and the strength of magnetic fields in intracluster 

■ environment and possibly will allow us detect the accretion shocks. 



1. Introduction 

Cosmic shocks emerge during structure formation in the universe, due to 
gravitationally driven supersonic gas infall onto collapsing objects. Astrophysical 
shocks are collisionless and capable of accelerating cosmic-rays (CRs) via first 
order Fermi mechanism [2]. CR electrons suffer severe energy losses dominated, 
in typical intra-cluster environment, by inverse Compton (IC) for energies above 
~ 150 MeV and Coulomb losses below that. Low energy, sub-relativistic CR 
protons also have short lifetimes due to Coulomb losses. However, for relativistic 
protons the dominant energy loss mechanism up to the energy threshold for photo- 
pair production is p-p inelastic collisions with a timescale longer than a Hubble 
time. Thus, once accelerated CR protons accumulate within large scale structures 
(LSSs) confined there by fiG strong turbulent magnetic fields [22]. 

It is of interest to investigate CR acceleration at LSS shocks for several 
reasons. Firstly, GCs exhibit non-thermal radiation. This mainly consists of 
diffuse radio emission extending on Mpc scales. It is thought to be synchrotron 



pp. 1-^ ©2002 by Universal Academy Press, Inc. 



2 



radiation but there is no consensus as to the origin of the emitting electrons. 
Also, excess of radiation with respect to thermal emission has been reported 
at hard X-ray [7] and, more controversially, extreme UV photon energies [10,4]. 
Secondly, if shocks acceleration operates efficiently, the proton component could 
bear a significant fraction of the total gas pressure, affecting the dynamics of 
cosmic structures [18]. Finally, CR electrons accelerated at inter-galactic shocks 
could contribute to the cosmic 7-ray background (CGB) [11]. This contribution 
addresses the role of 7-ray observations for investigating CR acceleration at LSS 
shocks. We describe the method in §^ the results in and conclude in §^ 

2. Numerical Simulation 

We a perform numerical simulation that models simultaneously the forma- 
tion of the LSS and the evolution of three CR species: namely shock accelerated 
protons and (primary) electrons and secondary e^ generated in p-p inelastic col- 
lisions of the simulated CR protons and the thermal gas. For the LSS we adopt 
a canonical, fiat ACDM model with total mass density Qm = 0.3, vacuum energy 
density = 0.7, Hubble constant Hq = 67 km s~^ Mpc"^, baryonic mass den- 
sity, Qb = 0.04, and spectral index and cluster normalization for initial density 
perturbations, Ug = 1 and ag = 0.9 respectively (see ref. [14] for full details). 

The CR dynamics, includes: injection at shocks, acceleration, energy losses 
and spatial transport. It is computed numerically through the code COSMOCR 
[13,14]. According to the adopted thermal leakage injection scheme, a fraction 
about 10~^ of the protons crossing the shock is injected as CRs. For shocks with 
Mach number, M = 4 — 10, responsible for most of the heating of the intergalactic 
gas [14], this roughly corresponds to converting 30% of the shock ram pressure into 
CR proton pressure. Further, we assume that the ratio of accelerated primary CR 
electrons to protons at relativistic energies is given by a parameter i?e/p ~ 10~^ 
[12,1]. In accord to the test particle limit of diffusive shock acceleration theory [2], 
the accelerated particles are distributed in momentum as a power-law with slope 
g = 4/(1 — M~^). All relevant loss mechanisms are accounted for [13,14]. We 
follow CR protons in the momentum range between 0.1 GeV/c and 10^ GeV/c, 
and CR electrons and e^ between 15 MeV/c and 20 TeV/c. 

3. Contribution to the Cosmic Gamma-ray Background 

Fig. |1] (left) shows the contribution to the CGB from CRs accelerated at 
LSS shocks [14]. It is computed as 




dz 



(1) 
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Fig. 1. Left: 7-ray background flux produced by cosmological CRs and EGRET data 
with tlieir error-bars (solid dots) [21]. Right: The IC 7-ray photon luminosity 
above 100 MeV from individual clusters as a function of the cluster X-ray core 
temperature, T^- EGRET experimental upper limits are from Reimer [19]. 



where e is the photon energy, j{e,z) is the computational-box-averaged spectral 
emissivity in units 'photons cm~ s at red-shift, attenuation factor 

due to photo-pair creation, 77 — >^ e^, c is the speed of light and Zmax an upper 
limit of integration. EGRET observational data (solid dots) [21] are also shown 
for comparison. We consider the following emission processes: IC emission of CR 
electrons scattering off cosmic microwave background photons (dot line), decay 
of neutral pions (tt*^ -f-f) produced in p-p inelastic collisions (dash line) and 
IC emission from secondary e^ (dot-dash line). In fig. |l] the total flux (solid 
line) corresponds roughly to a constant value at the level of 0.2 keV cm~^ s~^ 
sr~^ throughout the spectrum. It is dominated by IC emission from primary 
electrons. Fractions of order 30% and 10 % are produced by 7r°-decay and IC 
emission from secondary e^ respectively. 

All three components produce the same fiat spectrum, consistent with the 
observations. The result is a reflection of the fact that the CRs distributions 
were generated in strong shocks [14]. However, the computed flux is only ~ 15 
% of the observed CGB. It is difficult to imagine a higher contribution from vr^- 
decay and IC emission from e^ . In fact, if more CR protons were produced at 
shocks, CR-induced shock modifications would actually reduce the population of 
7-ray emitting protons (and e^ ). On the other hand, the fraction, rj, of shock 
ram pressure converted into CR electrons, can be constrained by comparing the 
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simulated clusters' 7-ray photon luminosity above 100 MeV to the upper limits 
set by the EGRET [20,19] for nearby GCs. This is done in fig. |1] (right panel). 
The simulation data are best-fit by the curve (solid line): 

L,(> 100 MeV) = 8.7 x 10« (j^) (^f ph s-. (2) 

According to our calculations [14], the EGRET upper limits require that t] < 
0.8%. This implies an upper limit on the computed 7-ray fiux of about 0.35 keV 
cm~^ s~^ sr~^ or a fraction of order ~ 25 % of the CGB. 

4. Synthetic emission maps and spectrum for a Coma-like prototype 

In this section we explore the spectral and spatial properties of 7-ray ra- 
diation between 10 keV and 10 TeV due to shock accelerated CRs in GCs. As 
non-thermal bremsstrahlung turns out unimportant [15], we consider the follow- 
ing emission processes: 7r°-decay and IC emission from both primary CR electrons 
and secondary e^ . We first compute the emission spectrum for the largest viri- 
alized object in the computational box which has an X-ray core temperature 
Tj. ~ 4 keV. The various emission components are then renormalized to the case 
of a Coma-like cluster. In particular, IC emission from electrons is rescaled ac- 
cording to equation The total number of e^ is rescaled assuming that these 
particles are responsible the for synchrotron emission of Coma radio halo 0. For 
the purpose we took a radio fiux at 1.4 GHz Siaghz = 640 mJy [5] and a two 
values for the volume average magnetic field (B), namely 1 and 3 fiG [9]. 

In the left quadrant of fig. |2] we present synthetic maps of the integrated 
photon fiux above 100 keV (top panels) and 100 MeV (bottom panels). They 
show that the emission from Tr'^-decay (bottom left) and e"*^ (top left) is confined 
to the cluster core where it creates a diffuse halo which rapidly fades with dis- 
tance from the center. In fact, e^ and ir^ are produced at the highest rate in 
the densest regions where both the parent CR ions and target nuclei are most 
numerous. On the other hand, because of severe energy losses, 7-ray emitting 
primary electrons are only found in the vicinity of strong shocks where they are 
accelerated. The emission (right panels) is extended and with a rich morphology 
refiecting the complex "web" of accretion shocks surrounding GCs [17]. The cor- 
responding synthetic spectra extracted from a core (top; with a 0.5° radius) and 
an outskirts region (bottom; a ring with inner and outer radii of 0.5° and 1.5° 
respectively) are shown in the central quadrant of fig. 0. They shows that the 
emitted radiation in the outskirts region is strongly dominated by IC emission 

*This assumption is not necessary but it is assumed for simplicity 
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Fig. 2. Left: Synthetic map of the integrated photon flux in units "ph cm~^ 
arcmin"^". Each panel measures 15 h~^Mpc on a side. Center: Synthetic spectra 
from 10 keV up to 10 TeV extracted from a core and an outskirts region. Right: 
Comparison of synthetic spectra with nominal sensitivity limits of future 7-ray 
observatories (thick-solid lines). Significance and observing time are respectively 
3c7 and 10*^ s for INTEGRAL-IBIS imagers (ISGRI and PICsIT), 5a and one year 
of all sky survey for EGRET and GLAST and 5a and 50 hour exposure for the 
Cherenkov telescopes. 



from primary primary e~. Conversely in the core region vr^'-decay (solid thin line) 
dominates at high photon energy (> 100 MeV) (top panel). Below ~ 100 MeV 
the relative contribution from primary and secondary e^ depends on the strength 
of the magnetic field. Finally, in the right quadrant the predicted radiation spec- 
tra in various spatial regions are compared to sensitivity limits^ of planned 7-ray 
observatories. 

5. Conclusions 

Clusters observations at 7-ray energy will provide important information 
about GCs. First, the extended emission is produced at the location of accretion 
shocks which, unlike weak merging shocks, have yet to be observed directly. Thus, 
detection and imaging of IC emission from primary electrons would provide an 
opportunity for their direct observation. The flux about 100 keV together with 
radio measurements, will allow an estimate of the intracluster magnetic field. It 
is important, however, to separate out the contribution from the CR electrons 
accelerated at the outlining shocks. In fact, because the magnetic field strength 
is expected to drop at the cluster outskirts, these electrons likely generate only a 
weak radio emission. Nonetheless, they produce a strong IC flux which can easily 



t These are for point sources and will be worse for extended sources considered here. 
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dominate the total soft 7-ray emission (cf. fig. Measuring the 7-ray flux at 
and above 100 MeV will allow us to confirm or rule out secondary e^ models for 
radio halos in GC [6,3,16] and to estimate the non-thermal CR pressure there 
[18]. In this perspective several authors estimated for nearby clusters the 7-ray 
flux from 7r°-decay. However, radiation flux from IC emission can be comparable 
[15]. Therefore, for a correct interpretation of the measurements the contributions 
from these two components will need to be separated as outlined here. 

Measuring the spectrum of IC 7-ray from primary CR electrons will also 
be instrumental for addressing their contribution to the CGB [11,14,8], the details 
of their acceleration mechanism and the physical conditions of the shocks. 
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